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1 Introduction 

The interactions of neutrinos with nucleons can provide valuable information on axial prop¬ 
erties and transition form factors [I]. For example, the nucleon’s axial form factor is rather 
badly known. It is usually reduced to a dipole ansatz, with one free parameter, the axial 
mass, remaining. This axial mass has been determined in many neutrino experiments on 
nucleons (or deuterons) and assumes a value of Ma ~ 1 GeV [2]. The dipole form, however, 
cannot really be constrained further by experiment [3] and indeed the vector form factors 
obtained from elastic electron scattering show a significantly more complicated dependence 
on the squared four-momentum [Ij. A similar situation exists for the transition form 
factors where even less is known. For example, for the A resonance the transition current 
involves 3 vector form factors and 3 axial ones. While the 3 vector form factors are reason¬ 
ably well determined by elektron-induced pion production on the nucleon, the 3 axial form 
factors are largely unknown. Present data [1] seem to be sensitive to only one of them and 
for it again the simple dipole form has been assumed. 

The investigation of interactions of neutrinos with nuclei may thus seem to be of only 
theoretical interest, since many elementary processes are not well understood. Indeed, 
there is a large theoretical interest in calculating the response of nuclei to external probes. 
This process can give valuable information about the nuclear many-body problem and the 
relevant reaction mechanisms, such as, e.g., the presence of interactions of the incoming 
probe with more than one nucleon UM- Such processes require theoretical methods that 
go beyond the so-called impulse approximation in which the interaction proceeds just with 
one nucleon at a time. It is encouraging to see that these investigations have made a huge 
progress over the last 20 years. First, the impulse approximation was extended to work not 
only with quasi-free nucleons, but instead with dressed nucleons that contain some effects 
of interactions with the surrounding nucleons in the nucleus |7]. These interactions can 
be partly summed up in ’spectral functions’ which are given by the imaginary part of the 
nucleon’s propagator inside the nucleus. The spectral functions contain effects of long-range 
mean-field potentials as well as effects of short-range correlations. Spectral functions have 
been very successfully employed in the description of the quasielastic response of nuclei 
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to electrons [8]. Not all of the many-body interactions can, however, be absorbed into a 
spectral function. The remaining interactions with neighboring nucleons then have to be 
treated explicitly by evaluating the many-body response through so-called 2 particle - 2 hole 
(2p2h) interactions [5l EJ [9l |T0l |TT] . Particularly noteworthy are recent ab-initio calculations 
of the quasielastic response of nuclei to incoming neutrinos that involves both vector and 
axial couplings m- Due to their huge computational requirements such calculations have 
become available only during the last few years; they are so far restricted to non-relativistic 
energies and to rather low energy transfers. At higher energies and energy transfers also 
inelastic excitations, either through nucleon resonance excitations or through deep inelastic 
scattering (DIS), take place and the treatment must be relativistic. 

While the investigation of interactions of neutrinos with nuclei may seem to be of 
only theoretical interest there is a very practical interest in such studies for long-baseline 
experiments that look for neutrino oscillations, such as, e.g., T2K, MINOS, NOvA and, 
in the future DUNE (formerly called LBNE). In such experiments the neutrino flux at a 
far detector is compared with that at a near detector. From that comparison the neutrino 
oscillation parameters, mixing angles and a possibly CP-invariance violating phase, can be 
extracted. What is actually compared is the event rate (flux times cross section at a given 
neutrino energy E^,) at a far detector with that at the near detector. The flux comparison 
thus requires the knowledge of the neutrino energy. The complication lies in the fact that 
the neutrino energy is not known because of the special production method of neutrinos 
as secondary decay products of particles, mostly pions and kaons, that were produced in 
primary reactions of protons with nuclei. The neutrino energy thus must be reconstructed 
event by event from the final state of the reaction. Two methods for this reconstruction are 
being considered: 

1. The first method is a so-called calorimetric method in which the energy of the hnal 
state particles is observed. If the detector were perfect this would give directly the 
incoming beam energy, through energy conservation. Real-life detectors, however, 
have limitations. They have acceptance thresholds and problems to observe certain 
particle classes, e.g. electrically neutral particles. The experiment then sees actually 
only a small part of the energy of the final state phase-space and must extrapolate 
from that to the full final state energy. For that extrapolation so-called neutrino 
generators have been used. Necessary theoretical input here is the knowledge of the 
initial neutrino-nucleon interaction and the hadron-hadron interactions in the hnal 
state. 

2. The second method is based on the fact that for quasielastic charged current scat¬ 
tering of a neutrino on a free neutron at rest the incoming neutrino energy can be 
determined completely from the outgoing lepton kinematics (energy, angle). This 
method obviously requires an experimental veto on any other particles - except for 
the proton - in the hnal state, for a correct identihcation of the reaction mechanism 
as being quasielastic scattering. The complication now comes because all ongoing and 
planned long-baseline experiments use nuclear targets (C, O, Ar, Fe), partly in order 
to increase the reaction rates and partly because of experimental safety considera¬ 
tions. In real life the neutrons are thus not free, but they are bound in a nucleus and 
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are Fermi-moving. Both of these facts lead to a smearing of the reconstructed energy 
around a sharp value, with an uncertainty width of about 60 MeV for a neutrino 
energy of about 1 GeV. 

While this presents a natural lower limit to the error with which the neutrino energy 
can be reconstructed, there is a more difficult problem to overcome when using nuclear 
targets that leads to significantly larger errors. This is the correct identification of the 
scattering event as being quasielastic. While for a free nucleon all that is needed is an 
experimental veto on outgoing hadrons (besides the proton), in a nuclear target this 
is not enough, because here now pions that were originally produced can be absorbed. 
The final state can then not be distinguished from one after QE scattering. Indeed, 
the misidentification of QE scattering has led to the extraction of unphysical values 
for the axial mass from experiments with nuclear targes names]. This also means 
that pion production, either through nucleon resonances or DIS, is always entangled 
with QE in a wide sense (true QE on one nucleon and 2p2h processes). In other 
words: QE cross sections can never be determined by purely experimental means, but 
always require an event generator to subtract the pion contributions. 

In order to get a feeling for the relative importance of QE scattering vs. pion production 
Eig.[T] shows these contributions for the flux expected at the near detector of the DUNE. 
It is noticeable that pion production processes (resonances and DIS) contribute about 2/3 
of the total cross sections and thus are the dominant component. Studies of the remainder, 
QE, then require first a quantitative understanding of pion production and an exact imple¬ 
mentation of that understanding in generators. This, in turn, requires knowledge both of 
the neutrino-nucleon pion production process, of in-medium changes to that cross section 
for bound nucleons and a very good description of pion-nucleon interactions. 

From these considerations it is clear that neutrino long-baseline experiments require 
a theoretical description of both the initial target ground state and the full hnal state; 
just inclusive cross sections are not enough. These exeriments thus need input from nu¬ 
clear structure and nuclear reaction theory. The former requires knowledge of the nuclear 
groundstate. The latter, on the other hand, requires knowledge of reaction mechanisms (1 
particle vs. 2 particle initial processes) as well as a state-of-the-art description of final state 
interactions. While for the former the knowledge of lepton-nucleus interactions is essen¬ 
tial, for the latter the essential requirement is knowledge of hadron-hadron interactions in 
the final state interactions. Without having these two areas quantitatively under control 
with state-of-the-art methods the goals of the precision era of neutrino physics can not be 
reached. 

2 Reaction mechanisms and oscillation signal 

We now have a closer look at the various reaction mechanisms contributing to the total 
neutrino-nucleus response in the LBNF beam: the nuclear target here is ^^Ar. All results 
in the following have been obtained within the transport theoretical GiBUU framework 
HZllIBI. GiBUU goes beyond Monte Carlo simulations in that it takes the nuclear potentials 
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Figure 1: (color online) Flux-averaged distribution {da/dQ'^) per nucleon for all events as 
a function of true for an target in the LBNF beam. The contribution labeled ’pions’ 
(solid blue line) gives the sum of all pion-producing processes (resonances, background and 
deep inelastic scattering), the one labeled ’QE’ (dashed red) depicts the sum of true one- 
body CCQE and of 2p2h processes. The solid black line gives the sum of both (taken 
from [T6]b 


into account. The code has been widely tested with the help of a broad range of nuclear 
reactions, with leptons, photons, hadrons and heavy ions as incoming beam. 

Fig. [2] shows the total event distribution and its decomposition into various production 
channels as a function of true neutrino energy. The true energies are those entering as 
input into the calculations. The overall energy-dependence of this event distribution is 
determined by the energy-distribution of the incoming neutrino beam. It is seen that - at 
the peak - true QE contributes about 1/3 of the total rate; the two next most important 
production channels are 2p-2h excitations and excitations of the A resonance. With a clear 
shift towards higher energies also DIS contributes on a similar level. The different processes 
thus exhibit different energy dependencies and they all overlap to a certain extent. This 
is a particular complication of experiments working with a neutrino flux in the few GeV 
region. At lower energies, e.g. the ones used by T2K, DIS and higher resonance excitations 
play no role so that essentially only A excitation contributes to pion production [20]. On 
the opposite, at significant higher beam energies larger than about 40 GeV the neutrino 
interactions are dominated by partonic degrees of freedom; there one has to deal with the 
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Figure 2: (Color online) Event distribution (normalized flux times cross section) per nucleon 
for LBNE vs. true energy (uppermost solid curve). Only events with 0 pions in the final 
state are taken into account. The various contributions to the total event rate are plotted 
as denoted in the figure (taken from |19jl. 


complications of in-medium effects on partons (EMC effect). This figure also illustrates the 
difficulties one faces in determining cross sections from such measured event rates which are 
products of flux and cross section. Any error made in reconstructing the neutrino energy 
will directly translate into an error of the cross section. 

In order to first get a feeling for the accuracy needed for the energy reconstruction in 
oscillation experiments a look at Eig. [^is helpful; this figure shows the expected oscillation 
signal for DUNE as a function as a function of neutrino Energy Ei, for some values of of two 
neutrino properties, mixing angle 613 and the CP-violating phase 5cp- The three curves 
under the flux profile can be distinguished from each other only if the neutrino energy can 
be determined to better than about 100 MeV. This gives a first hint at the accuracy that is 
needed at DUNE. Errors in the energy reconstruction due to event-misidentification cause 
not just a shift of the energy axis, but instead distort the whole event distribution [2^[23] . 

For the neutrino oscillation analysis the neutrino energy plays an essential role because 
it enters into the oscillation formula. Errors in the reconstructed energy thus have an impact 
on error in the mixing angles and phases. This is illustrated in Fig. which in its lower two 
curves shows the expected event distribution for electron appearance at the far detector 
in the LBNF beam. Here, the oscillation signal is plotted for two different event samples 
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Figure 3: Ve appearance probability at a distance of 1300 km calculated for standard os¬ 
cillation mixing angles. The four colored curves illustrate the sensitivity of the expected 
signal to the neutrino mixing angle @13 and the CP-violating phase 5cp- The black peak 
shows the expected energy distribution for the neutrino beam (taken from j 21 ]i. 


as a function of true neutrino energy (solid curves) as well as, on the other hand, of an 
energy reconstructed from the outgoing electron kinematics assuming a true QE process 
(dashed curves). The upper two curves give the oscillation signal that is obtained from an 
event sample with 0 pions. These two curves, calculated with GiBUU, directly correspond 
to the red curve in Fig. The event rate vs. reconstructed energy (upper dashed curve) 
is distorted as compared to the one vs. true energy (upper solid curve) and shifted by 
more than 500 MeV in its maximum even though only events with 0 outgoing pions have 
been used for the reconstruction. This is clearly above the accuracy required to distinguish 
between the various parameter scenarios in Fig. 

A drastic improvement happens when the event sample is further restricted to contain 
one and only one proton (plus any number of neutrons). Now the difference between 
the lower two curves is at most 100 MeV. Since the energy reconstruction is based on 
the dynamics of a true ( 1 -body) process this implies that requiring 1 proton in addition 
to 0 pions gives a cleaner identihcation of true QE. This result depends crucially on the 
theoretical description of final state interactions, not just between pions and nuclei, but also 
between nucleons. The final state cascade must allow for an avalanche effect in which in most 
cases the outgoing proton from a neutrino-induced charged current scattering rescatters thus 
increasing the multiplicity of final state nucleons and and simultaneously lowering their 
energies. An event with only 1 hnal state proton then is a rather clean, direct signal of QE 
which has not undergone this final state scattering. This has recently also been exploited 
in a study of QE scattering by the MINERi^A experiment |24j . 
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Figure 4: Vf. appearance event distribntion (normalized flnx times cross section) per nncleon 
for LBNE vs. true (solid curve) and reconstrncted (dashed cnrve) energy. The upper two 
curves show the results obtained from an event sample with 0 pions, the two lower curves 
are obtained from a sample with 0 pions, 1 proton and X neutrons (taken from [H]). 


The same improvement also shows np for the difference between the true and the recon¬ 
structed oscillation signal in dependence on bcp m. Experiments looking for this phase 
would be well advised to look at events with 0 pions, 1 proton and X (unobserved) nentrons. 

3 Summary 

Neutrino oscillation parameters can be extracted only if the incoming nentrino energy is 
known. The latter has to be reconstructed from final state particles. This reconstruction 
requires knowledge of nentrino-nucleon interaction rates in medium and of final state inter¬ 
actions between the outgoing hadrons. It also requires a theory that is able to describe the 
complete time-development of the nentrino-nucleus reaction; just inclusive cross sections 
are not enough. Employing state-of-the-art theoretical methods of nuclear physics, both for 
the description of the ground state and of the reaction mechanisms, is essential to be able 
to take event generators out into new regions of energy and target mass. 

Over the last few years tremendous progress has been made in this respect. The theoreti¬ 
cal methods are available and have also been ~ at least partly - implemented into theoretical 
descriptions of a neutrino-nuclens reaction. It has also been realized that presently used 
event generators often lag behind in their implementation of present-day’s nnclear physics. 
However, it is clearly disconcerting for any nuclear theorist to see that even in very recent 
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experiment analyses by leading experiments still outdated methods with unphysical param¬ 
eters are being used. One of the two most glaring examples are the large, unphysical axial 
mas of Ma ~ 1.25 GeV that still underlies the description of QE scattering in recent exper¬ 
imental analyses as well as the still ongoing use of the Rein-Sehgal formfactors for nucleon 
resonances that are known to fail in their description of electroproduction data [25t I26j . 
Both of these severe shortcomings have been realized now for about 6 years, but still have 
not caused any modification of the generators used by these experiments. This is a clearly 
very unsatisfactory situation, both from a nuclear theory point of view and from a neutrino 
experimental point of view, since it runs counter to the ambitions of the precision era of 
neutrino physics that has now begun. How uncertainties in the generators used actually 
affect the oscillation mixing angles has been illustrated in recent work |27l [28] . 
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